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Acronyms: 

HRS  hepatocyte  growth  factor-regulated  tyrosine  kinase  substrate 

NF2  Neurofibromatosis  2 

RT4  cells  Schwannoma-like  cells  obtained  by  ENU  mutagenesis 

EGF  epidermal  growth  factor 

ST AT  Proteins  Signal  transducers  and  activators  of  transcription 
MEF  Mouse  embryonic  fibroblasts 

Introduction: 

We  would  like  to  make  it  clear  at  the  outset  that  we  are  more  than  willing  to  deposit 
our  raw  microarray  data  into  a  database  should  the  Army  Medical  Branch  have  any 
suggestions  for  such  a  site  maintained  by  the  army.  These  data  may  be  of  interest  to 
other  investigators  in  the  oncology  field. 

The  NF2  protein,  also  called  schwannomin  or  merlin,  is  a  member  of  the  ezrin-radixin- 
moesin  family  of  proteins.  Based  on  mutation  analysis  and  protein  expression  studies,  it  is 
thought  that  NF2  is  a  tumor  suppressor  gene  and  that  biallelic  inactivation  is  required  for 
phenotypic  expression.  Studies  in  cells  that  are  NF2  deficient  or  that  overexpress  NF2 
cDNAs  have  shown  a  powerful  role  of  NF2  in  regulation  of  cell  morphology  and  proliferation. 
Despite  these  advances  little  is  known  about  the  action  of  NF2  in  specific  signaling 
pathways.  This  knowledge  is  of  importance  because  it  may  lead  the  way  to  identifying 
treatments  that  modify  these  pathways. 

Expression  profiling  is  a  powerful  novel  technique  to  examine  changes  in  the  expression  of 
a  large  number  of  genes  at  the  same  time.  Parallel  analysis  of  gene  expression  reflects  the 
changing  view  of  signaling  pathways  towards  signaling  networks  that  have  multiple  and 
complex  feedback  and  feed  forward  loops.  Different  phenotypic  states  of  a  cell  can  be 
translated  into  specific  gene  expression  signatures.  Expression  profiling  is  not  restricted  to 
known  genes,  but  changes  in  the  expression  of  genes  without  known  function  can  also  be 
detected. 

We  are  proposing  to  use  gene  expression  profiling  to  determine  changes  in  gene  expression 
following  regulated  expression  of  NF2  cDNAs  in  schwannoma  and  meningioma  cells.  We  will 
test  the  following  hypotheses:  1.  Regulated  expression  of  NF2  results  in  specific  changes  in 
the  expression  profile  at  specific  time  points.  2.  A  subset  of  regulated  transcripts  differs 
between  cells  expressing  NF2  isoform  1  and  NF2  isoform  2.  3.  NF2  induction  results  in 
partially  shared  expression  changes  in  different  cell  types.  4.  Expression  changes  are 
context  dependent;  expression  profiles  in  confluent  cells  are  different  from  non-confluent 
cells. 


Body 

Statement  of  Work  Year  1 

In  the  first  year  we  will  begin  with  the  analysis  of  timed  NF2  expression.  In  the  second  half  of  the 
first  year  we  will  examine  the  effects  of  overexpressing  NF2  iso  form  1  or  isoform  2  in  rat  RT4 
schwannoma  cells. 
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For  those  genes  that  show  consistent  changes  in  expression  we  will  perform  quantitative  northern 
blot  analysis  to  validate  expression  changes. 

We  will  continue  the  generation  of  meningioma  cells.  Lines  expressing  wiidtype  and  mutant  forms  of 
NF2  will  be  generated  in  the  Tet/On  system. 

We  have  established  parameters  for  the  induction  of  NF2  expression  in  RT4  cells.  We  had 
problems  generating  meningioma  cells  that  stably  express  NF2.  We  therefore  used  mouse 
embryonic  fibroblasts  (MEF).  cDNAs  encoding  NF2  isoforms  1  and  2  were  put  under  the 
control  of  the  tet-responsive  promoter  and  introduced  into  MEF/3T3  TET-OFF  cells  using 
retroviral  transduction. 

These  cells  show  a  robust  induction  of  NF2  expression  at  the  RNA  and  protein  levels  after 
removal  of  doxycycline. 

We  have  purified  RNAs  from  RT4  and  MEF  cells  and  used  these  for  expression  profiling  with 
the  Affymetrix  MG_U74Av2  chips.  At  a  significance  level  of  0.05,  we  found  that 
approximately  50  genes  were  regulated  upon  NF2  expression.  This  was  true  for  either 
isoform.  Interestingly,  no  regulated  genes  were  shared  between  isoform  1  and  2. 

We  have  now  completed  this  study  with  the  final  conclusion  that  the  overall  number  of 
regulated  genes  was  less  than  the  number  expected  by  chance.  Thus,  we  concluded  that  in 
the  given  experimental  paradigm,  expression  of  NF2  did  not  induce  any  measurable  effects 
on  gene  expression.  These  results  are  presented  in  detail  in  the  appended  manuscript  that 
is  now  published  in  the  Journal  of  Cellular  Biochemistry. 


Statement  of  work  year  2 
Year  2 

We  will  evaluate  gene  expression  profiles  in  meningioma  cells.  Unless  software  updates  are 
implemented  at  that  time,  we  expect  that  comparison  of  human  and  rat  profiles  will  be  time- 
consuming  and  may  have  to  be  done  manually. 

For  those  genes  that  show  consistent  expression  changes  in  meningioma  cells  we  will  perform 
quantitative  northern  blot  analysis  to  validate  expression  changes. 

Based  on  the  results  shown  in  our  manuscript,  we  did  not  think  that  it  was  warranted  to  apply  the 
same  experimental  paradigm  to  the  study  of  meningeal  cells.  We  have  therefore  changed  our 
approach  of  detecting  transcriptional  changes  making  use  of  our  recent  findings  that  expression  of 
NF2  and  its  binding  protein  HRS  influence  STAT-signaling  in  response  to  EGF  (epidermal  growth 
factor).  STATS  (Signal  transducers  and  activators  of  transcription)  are  proteins  that  are  important  in 
intracellular  signaling.  STATs  3  and  5  have  recently  been  linked  to  EGF  signaling  in  tumor  cells. 

With  the  remaining  funds  we  have  now  analyzed  the  effect  of  EGF  treatment  in  RT4  cells  with  and 
without  induction  of  Hrs.  Quite  surprisingly,  relatively  few  studies  have  examined  the  effects  of 
induction  of  EGF  treatment  on  cells  at  the  microarray  level.  These  studies  also  gave  us  information 
on  the  changes  in  gene  expression  induced  by  long-term  overexpression  of  Hrs. 

The  experiment  was  carried  out  in  the  following  way: 

RT4  cell  lines,  which  express  HRS  upon  addition  of  doxycycline  (Tet-On),  were  incubated  in  DMEM 
containing  10  %  Tet-approved  FBS  (Clontech),  2  yg/ml  doxycycline,  100  pg/ml  G418,  and  250 
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|jg/ml  hygromycin  B  until  they  reached  50%  confluence.  Those  cell  lines  were  washed  three  times 
with  DMEM  and  incubated  again  in  5ml  DMEM. 

For  one  set  of  cell  lines,  lpg/ml  doxycycline  was  added  to  induce  Hrs  expression,  the  other  set  was 
treated  with  vehicle.  Then,  lOOmg/ml  of  EGF  was  added  to  both  sets  of  cell  lines.  The  cell  lines  were 
harvested  at  lOmin,  30min,  1  hour,  and  2  hours  after  EGF  treatment  and  total  RNAs  were  purified. 
The  cell  lines  before  EGF  treatment  were  also  harvested  as  reference.  The  RNAs  were  subjected  to 
the  microarray  experiments  using  the  RU230A  GeneChip  (Affymetrix).  Array  images  (CEL  files  from 
the  GENECHIP  program,  Affymetrix)  were  imported  into  the  DCHIP  program  and  the  expression 
levels  were  calculated  using  the  model-based  method  after  normalization  with  DCHIP's  rank 
invariant  protocol. 


Hrs  ♦ 

Hrs  - 

Fig.  1:  Expression  changes  for  the  Early  Growth 
Response-1  transcript,  one  of  the  10  most 
activated  genes  by  EGF  treatment  in  RT4  cells. 
Early  Growth  Response- 1  transcript  was  typical  for 
a  number  of  transcripts  showing  maximum 
induction  reached  at  30  min.  This  gene  also 
showed  a  significant  effect  of  reduced  expression 
at  30  min  in  the  presence  of  Hrs. 


All  transcripts  studied  showed  an  increase  in  expression.  Several  genes  showed  expression  changes 
in  the  order  of  10  to  15  fold  upon  EGF  induction,  most  of  them  within  lOmin  and  maintained 
through  60  min. 


Hrs  effect: 

It  was  also  interesting  to  compare  the  effect  of  Hrs  expression  in  cells  treated  with  EGF  or  with 
vehicle.  HRS  was  overexpressed  in  serum  starved  condition,  and  the  RNAs  were  taken  18  hours 
after  dox-treatment.  Table  1  list  the  genes  that  showed  significant  or  highly  significant  changes  in 
expression  when  RT4  cells  expressing  high  steady-state  levels  of  Hrs  were  compared  with  cells 
expressing  only  endogenous  Hrs  at  relatively  low  levels.  This  comparison  was  performed  in  the 
presence  and  in  the  absence  of  EGF. 

Only  four  genes  showed  reduced  expression  with  HRS  expression.  One  of  these  genes, 
glutaredoxin-1,  was  represented  by  two  independent  probe  sets.  The  two  probe  sets  gave  almost 
identical  results  attesting  to  the  quality  of  our  microarray  experiments. 

Interestingly,  for  most  genes,  whose  expression  levels  were  altered  by  Hrs,  the  relative  reduction 
was  maintained  in  the  presence  or  absence  of  EGF.  The  genes  fall  into  several  classes  including 
apoptosis,  metabolic,  t-RNA  ligases,  and  cytokines/growth  factors.  Of  note,  one  of  the  genes, 
syndecan,  is  known  to  also  interact  with  the  known  NF2  binding  protein  syntenin. 
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Table  1:  Long-term  effects  of  Hrs  expression  on  gene  expression  in  the  presence  or  absence  of  EGF 


Gene 

Protein  Function 

EGF  treated 

EGF  untreated 

Fold  change 

P  value 

Fold  change 

P  value 

glutaredoxin  1  (thioltransferase) 

Metabolic 

-1.8 

0.0014 

-1.57 

0.0214 

neurofilament,  light  polypeptide 

Structural 

-1.7 

0.0032 

-1.45 

0.0492 

neurofilament,  light  polypeptide 

Structural 

-1.65 

0.0033 

-1.45 

0.0121 

glutaredoxin  1  (thioltransferase) 

Metabolic 

-1.59 

0.0018 

-1.41 

0.0318 

ESTs,  Highly  similar  to  CAV1_MOUSE  Caveolin-1 
[M.musculus] 

Endocytosis 

-1.52 

0.0121 

-1.41 

0.0380 

ESTs,  Highly  similar  to  A33267 
methylenetetrahydrofolate  dehydrogenase  (NAD+) 

Metabolic 

1.42 

0.0008 

1.55 

0.0016 

l-gicerin 

Receptor 

1.44 

0.0362 

1.69 

0.0160 

interferon-related  developmental  regulator  1 

GF/Cytokine 

1.47 

0.0316 

1.69 

0.0030 

ESTs,  Moderately  similar  to  JC5224  methionine--tRNA 
ligase  (EC  6.1.1.10)  -  human  [H.sapiens] 

tRNA  ligase 

1.48 

0.0098 

1.55 

0.0039 

Syndecan 

Receptor;  Binds  the 
NF2  binder  syntenin 

1.49 

0.0262 

1.59 

0.0049 

ESTs,  Highly  similar  to  growth  hormone  inducible 
transmembrane  protein;  EST  C77840  [Mus  musculus] 

Receptor 

1.52 

0.0039 

1.41 

0.0048 

ESTs,  Weakly  similar  to  SYD  RAT  ASPARTYL-TRNA 
SYNTHETASE  (ASPARTATE-TRNA  LIGASE) 

Metabolic 

1.54 

0.0018 

1.74 

0.0159 

EST  unknown 

1.55 

0.0139 

1.67 

0.0156 

ESTs,  Highly  similar  to  S50053  tryptophan-tRNA  ligase 
(EC  6. 1.1. 2)  alpha-2  chain  -  mouse  [M.musculus] 

tRNA  ligase 

1.55 

0.0003 

1.58 

0.0005 

putative  fatty  acid  desaturase 

Metabolic 

1.55 

0.0090 

1.56 

0.0030 

ESTs,  Highly  similar  to  SYAJHUMAN  Alanyl-tRNA 
synthetase  (Alanine--tRNA  ligase)  (AlaRS)  [H.sapiens] 

tRNA  ligase 

1.57 

0.0020 

1.43 

0.0129 

ESTs,  Highly  similar  to  SYC  HUMAN  CYSTEINYL- 
TRNA  SYNTHETASE  (CYSTEINE-TRNA  LIGASE) 
(CYSRS)  [H.sapiens] 

tRNA  ligase 

1.58 

0.0049 

1.65 

0.0014 

sequestosome  1 

proteasome 

1.59 

0.0035 

1.76 

0.0014 

G  protein-coupled  receptor  56 

Receptor 

1.6 

0.0263 

1.51 

0.0265 

ESTs,  Highly  similar  to  A55314  glycine-tRNA  ligase 
(EC  6.1.1.14)  precursor  [validated]  -  human  [H.sapiens] 

tRNA  ligase 

1.63 

0.0008 

1.7 

0.0052 

ESTs,  Weakly  similar  to  GLYC_MOUSE  Serine 
hydroxymethyltransferase 

Metabolic 

1.63 

0.0005 

1.79 

0.0002 

ESTs,  Highly  similar  to  type  1  transmembrane  protein 
Fn14;  fibroblast  growth  factor  regulated  protein  2  [Mus 
musculus] 

Receptor 

1.64 

0.0126 

1.62 

0.0181 

gb:AB020967.1  /DB XREF=gi:4827158  /GEN=NIPK 

Apoptosis 

1.69 

0.0013 

1.71 

0.0006 

activating  transcription  factor  5 

TF 

1.7 

0.0008 

1.84 

0.0329 

ESTs,  Moderately  similar  to  SYEPJHUMAN  Bifunctional 
aminoacyl-tRNA  synthetase 

Metabolic 

1.73 

0.0038 

1.65 

0.0332 

ubiquitin  C 

proteasome 

1.74 

0.0145 

1.9 

0.0033 

activating  transcription  factor  ATF-4 

TF 

1.81 

0.0035 

1.96 

0.0003 

growth  differentiation  factor  15 

GF/Cytokine 

1.87 

0.0274 

1.78 

0.0203 

nuclear  protein  1 

Apoptosis 

1.87 

0.0004 

1.8 

0.0011 

gb:AB020967.1  /DB XREF=gi:4827158  /GEN=NIPK 

Apoptosis 

1.94 

0.0004 

1.81 

0.0002 

Asparagine  synthetase 

Metabolic 

2.22 

0.0000 

2.12 

0.0010 

CCAAT/enhancer  binding  protein  (C/EBP),  beta 

TF 

2.41 

0.0010 

1.94 

0.0321 

ESTs,  Weakly  similar  to  interferon,  alpha-inducible 
protein  27-like  [Rattus  norvegicus]  [R.norvegicus] 

GF/Cytokine 

3.44 

0.0000 

3.2 

0.0201 
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Key  Research  Accomplishments 

•  Establishment  of  cell  lines  that  express  NF2  under  control  of  the  TET-regulator. 

•  Establishment  of  array  profiles  for  both  NF2  isoforms. 

•  Establishment  of  a  paradigm  of  STAT-activation  in  response  to  EGF. 

•  Gene  expression  profiles  of  Schwannoma-like  RT4  cells  under  conditions  of  NF2 
expression,  Hrs  expression,  and  treatment  with  EGF. 


Reportable  Outcomes 

Manuscripts  published  : 

Presentations 

MK  Oh  &  Pulst  "The  analysis  of  oligonucleotide  microarray  data" 

DNA  Microarray  Workshop  II,  University  of  California,  Los  Angeles,  CA.  USA  Sep.  16-17,  2002 
MK  Oh  "DNA  microarray:  A  promising  tool  for  pathway  engineering" 

Department  of  Chemical  and  Biological  Engineering,  Korea  University.  South  Korea,  May.  29,  2003 

Personnel  supported  by  the  award 

Oh  MK  Post-doctoral  Fellow 

Pulst  SM  Principal  Investigator 


Conclusions 

Gene  expression  profiling  is  still  in  its  infancy.  We  are  beginning  to  answer  some  basic 
questions  regarding  the  use  of  stably  transfected  cell  lines  and  the  role  of  schwannomin  and 
HRS.  We  have  now  established  that  quite  to  our  surprise  changes  in  steady-state  levels  of 
NF2  or  HRS  show  relatively  few  effects  on  the  expression  of  many  genes. 

We  altered  our  strategy  and  have  examined  changes  in  the  expression  profile  acutely  after 
treatment  with  growth  factors.  These  results  will  provide  a  framework  to  evaluate  the 
alterations  seen  in  human  tumors  associated  with  NF2  mutation. 
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Genetic  Heterogeneity  of  Stably  Transfected  Cell 
Lines  Revealed  by  Expression  Profiling  With 
Oligonucleotide  Microarrays 

Min-Kyu  Oh/  Daniel  R.  Scoles,1  Carrie  Haipek,2  Andrew  D.  Strand/  David  H.  Gutmann/2 
James  M.  Olson/  and  Stefan-M.  Pulst1'4* 

’Rose  Moss  Ncurogenetics  Laboratory  for  Parkinson  and  Related  Diseases,  CSMS  Burns 
and  Allen  Research  Institute,  Division  of  Neurology,  Cedars-Sinai  Medical  Center, 

8700  Beverly  Boulevard,  Los  Angeles,  California  90048 

department  of  Neurology,  Washington  University  of  Medicine,  St.  Louis,  Missouri  63110 
3Fred  Hutchison  Cancer  Center,  Seattle,  Washington  98109 

department  of  Medicine  and  Neurobiology,  David  Ceffen  School  of  Medicine  at  UCLA, 
Los  Angeles,  California 


Abstract  Large-scale  gene  expression  measurements  with  oligonucleotide  microarrays  have  contributed 
tremendously  to  biological  research.  However,  to  distinguish  between  relevant  expression  changes  and  falsely  identified 
positives,  the  source  and  magnitude  of  errors  must  be  understood.  Here,  we  report  a  source  of  biological  variability  in 
microarray  experiments  with  stably  transfected  cell  lines.  Mouse  embryonic  fibroblast  (MEF/3T3)  and  rat  schwannoma 
(RT4)  cell  lines  were  generated  to  provide  regulatable  schwannomin  expression.  The  expression  levels  of  29  samples  from 
five  different  mouse  embryonic  fibroblast  clonal  cell  lines  and  18  samples  from  3  RT4  cell  lines  were  monitored  with 
oligonucleotide  microarrays.  Using  hierarchical  clustering,  we  determined  that  the  changes  in  gene  expression  induced 
by  schwannomin  overexpression  were  subtle  when  compared  with  those  detected  as  a  consequence  of  clonal  selection 
during  generation  of  the  cell  lines.  The  hierarchical  clustering  implies  that  significant  alterations  of  gene  expression  were 
introduced  during  the  transfection  and  selection  processes.  A  total  of  28  genes  were  identified  by  Kruskal -Wall is  rank  test 
that  showed  significant  variation  between  clonal  lines.  Most  of  them  were  related  to  cytoskeletal  function  and  signaling 
pathways.  Based  on  these  analyses,  we  recommend  that  replications  of  experiments  with  several  selected  cell  lines  are 
necessary  to  assess  biological  effects  of  induced  gene  expression.  J.  Cell.  Biochem.  90:  1068-1078,  2003. 
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Expression  profiling  with  oligonucleotide 
microarrays  has  greatly  influenced  biological 
research.  By  measuring  expression  levels  of 
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thousands  of  genes  simultaneously,  microarray 
experiments  facilitate  the  study  of  complicated 
phenotypes,  such  as  cancer  classification  or 
circadian  biological  clocks  [Golub  et  al.,  1999; 
Storch  et  al.,  2002 1.  This  technology  has  also 
been  employed  to  identify  genes  in  cell  culture 
model  systems  that  reflect  the  function  of 
transcriptional  regulators,  such  as  p53  [Zhao 
et  al.,  2000],  EGR1  [Svaren  et  al.,  2000],  E2F 
[Ma  et  al.,  2002],  and  Pax3  [Mayanil  et  al., 
2001],  or  to  characterize  the  effects  of  patho¬ 
genic  protein  expression,  such  as  BRCA1 
IWelcsh  et  al.,  2002]  and  huntingtin  [Sipione 
et  al.,  2002]. 

Because  the  number  of  genes  simultaneously 
studied  is  large  in  microarray  experiments,  it 
is  important  to  characterize  the  sources  and 
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magnitudes  of  error  in  an  effort  to  sort  false 
positives  and  negatives  from  truly  regulated 
genes.  Systematic  analyses  of  errors  inherent  in 
oligonucleotide  microarray  experiments  have 
just  begun.  Generally,  ‘technical  (or  experi¬ 
mental)  variability’  is  defined  as  the  errors 
incorporated  during  experimental  steps,  such 
as  chip-to-chip  variations,  mRNA  purification, 
cDNA  probes  synthesis,  hybridization,  scan¬ 
ning,  and  image  analysis.  Although  this  techni¬ 
cal  variability  is  reduced  by  the  standardization 
of  experimental  protocols,  a  substantial  level  of 
noise  can  still  occur  during  experimental  pro¬ 
cedures.  A  recent  study  with  oligonucleotide 
microarray  experiments  showed  that  the  hybri¬ 
dization  step  is  the  major  source  of  the  technical 
variability.  The  hybridization  noise  is  highly 
dependent  on  the  expression  level  of  the  genes, 
suggesting  that  the  significance  of  gene  expres¬ 
sion  fold  changes  must  be  adjusted  depending 
on  the  gene  expression  levels  |Tu  et  al.,  2002]. 

In  contrast,  ‘biological  variability’  originates 
from  different  sources  of  RNA,  such  as  differ¬ 
ent  tissue  samples  in  a  replicated  experiment. 
In  general,  the  biological  variability  is  signifi¬ 
cantly  higher  than  the  technical  variability 
[Bakay  et  al.,  2002;  Novak  et  al.,  2002 1.  For 
example,  significant  background  variations  in 
expression  levels  were  observed  among  mice 
that  were  genetically  identical  and  had  been 
housed  under  the  same  conditions.  Many  of 
these  differentially  expressed  genes  in  those 
mice  were  heat  shock,  immune  responsive,  and 
hormone-regulated  genes  [Pritchard  et  al., 
2001;  Novak  et  al.,  2002].  In  other  experiments, 
the  expression  levels  varied  significantly  when 
the  samples  were  taken  from  different  regions 
of  the  same  muscle  tissue  | Bakay  et  al., 
2002].  These  results  suggest  that  the  biological 
samples  even  in  replicated  experiments  may 
contain  significant  background  variation  in 
gene  expression  levels. 

Here,  we  report  significant  biological  varia¬ 
bility  detected  in  stably  transfected  clonal  cell¬ 
lines.  With  oligonucleotide  microarrays,  we 
monitored  the  expression  level  changes  of 
several  different  clones  selected  for  the  regu¬ 
lated  expression  of  the  neurofibromatosis  2 
(NF2)  tumor  suppressor  protein,  schwannomin 
(or  merlin).  Hierarchical  clustering  of  the  ex¬ 
pression  data  showed  that  the  expression 
profiles  from  different  clonal  lines  were  signi¬ 
ficantly  different  and  these  differences  were  not 
significantly  influenced  by  expression  of  the 


transgene.  Our  results  suggest  that  significant 
gene  expression  heterogeneity  is  introduced 
during  the  process  of  selecting  stably  trans¬ 
fected  clones.  This  variation  may  be  greater 
than  that  resulting  from  expression  of  the 
transgene  of  interest. 

MATERIALS  AND  METHODS 
Cells  and  Transfection 

Full-length  cDNAs  coding  for  schwannomin 
isoforms  1  and  2 1  Scoles  et  al.,  1998]  were  cloned 
into  the  pRevTRE  plasmid  (Clontech)  under 
the  Tet-inducible  promoter.  The  plasmids  were 
transfected  into  a  packaging  cell  line  PT67 
(Clontech)  using  Superfect  reagent  (Qiagen). 
The  retrovirus  produced  from  the  packaging  cell 
lines  was  infected  into  mouse  embryonic  fibro¬ 
blast  Tet-Off  Cell  line  (Clontech).  Hygromycin 
resistant  Tet-Off  MEF/3T3  clones  were 
selected,  diluted,  and  plated  in  96-well  plates 
to  isolate  single  cell  clones.  The  cells  were  in¬ 
cubated  until  confluent  and  seeded  into  three 
96-well  plates,  two  plates  for  testing  schwanno¬ 
min  inducibility  by  dox/no  dox  comparison  and  a 
third  for  propagation.  From  200  colonies,  four 
cell  lines  (15-7  and  4-6  for  isoform  1,  and  23-6 
and  2-19  for  isoform  2)  showing  the  best  induc¬ 
tion  and  suppression  of  schwannomin  expres¬ 
sion  in  dox/no  dox  comparisons  were  chosen  for 
further  experiments.  The  four  selected  cell  lines 
and  the  parental  MEF/3T3  Tet-Off  cell  line  were 
maintained  in  DMEM  containing  10%  Tet- 
approved  FBS  (Clontech),  2  pg/ml  doxycycline, 
100  pg/ml  G418,  and  250  pg/ml  Hygromycin  B. 

The  construction  and  maintenance  of  induci¬ 
ble  Tet-On  RT4  cell  lines  for  schwannomin  and 
missense  mutant  (L64P)  schwannomin  were 
described  previously  [Gutmann  et  al.,  2001]. 

RNA  Preparation  and  cDNA  Probe  Synthesis 

MEF/3T3  cell  lines  were  incubated  at  50% 
confluence  in  100-mm  plate  and  washed  twice 
with  37°C  dox-free  media.  Cells  were  reincu¬ 
bated  for  24  h  in  medium  with  or  without 
doxycycline  for  repression  and  induction  of 
schwannomin  expression.  Total  RNAs  were 
purified  from  the  cells  using  the  RNeasy  Mini 
Kit  (Qiagen)  according  to  the  manufacturer’s 
protocol.  For  immunoblotting  experiments, 
10%  of  cell  lysates  were  saved.  Labeled  cRNA 
probes  were  generated  from  10  pg  of  total  RNA 
and  subjected  to  microarray  experiments  using 
the  MG_U74Av2  GeneChip  (Affymetrix)  as 
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described  previously  [Luthi-Carter  et  al.,  2002]. 
For  minimizing  variability  during  the  experi¬ 
ments,  all  cell  lines  were  grown  under  the  same 
conditions  at  the  same  time.  The  RNA  purifica¬ 
tion  and  chip  experiment  were  performed  at  the 
same  time  by  the  same  person. 

RT4  cell  lines  were  incubated  in  media 
containing  tetracycline-free  fetal  bovine  serum 
(Clontech)  until  they  reached  50%  confluence 
and  1  pg/ml  doxycycline  was  added  into  the 
media  to  induce  schwannomin  expression.  Total 
RNAs  were  purified  at  0,  6,  and  24  h  after  the 
addition  of  doxycycline  and  subjected  to  the 
microarray  experiments  using  RG  U34A  Gen- 
eChip  (Affymetrix). 

Data  Analysis 

Array  images  (CEL  files  from  the  GENECHIP 
program,  Affymetrix)  were  imported  into  the 
DCHIP  program  and  the  expression  levels  were 
calculated  using  the  model-based  method  after 
normalization  with  DCHIP’s  rank  invariant 
protocol  [Li  and  Wong,  2001].  Gene  filtering 
based  on  the  intensities  or  ‘present’  call  and 
Student’s  /-test  were  performed  using  the 
DCHIP  program.  The  /-test  was  also  performed 
with  the  expression  levels  calculated  by  the 
Affymetrix  GENECHIP  program.  The  results 
were  very  similar  to  those  obtained  using  the 
DCHIP  protocol  (data  not  shown). 

For  cluster  analysis  and  Kruskal- Walhs 
rank  test,  the  expression  levels  from  DCHIP 
were  divided  by  the  average  expression  level  of 
the  respective  gene  in  the  entire  set  of  chips 
(26  chips  for  MEF/3T3  cell  line  and  17  chips 
for  RT4  cell  line).  The  expression  ratios  were 
imported  into  the  CLUSTER  program  lEisen 
et  al.,  19981,  log- transformed  and  analyzed  by 
the  average  linkage  clustering  method.  The 
clustered  tree  was  drawn  by  the  TREEVIEW 
program  [Eisen  et  al.,  1998].  The  log- transform¬ 
ed  expression  ratios  were  used  for  Kruskal- 
Wallis  rank  test,  a  non-parametric  analysis 
method  similar  to  one-way  analysis  of  variance 
(ANOVA)  to  reveal  the  differentially  expressed 
genes  in  different  cell  line  clones.  Kruskal- 
Wallis  rank  test  was  conducted  using  the  S- 
Plus2000  program  (Insightful  Corp.). 

Immunoblotting 

The  cell  lysates  from  MEF/3T3  cell  line  sam¬ 
ples  were  precipitated  by  adding  five  volumes 
of  acetone  (-20°C),  incubating  for  15  min  at 
-20°C,  and  centrifuging  for  20  min  at  4°C.  The 


pellets  were  air-dried  and  dissolved  in  20  pi 
sample  buffer  consisting  of  50  mM  Tris-HCl 
(pH  6.8),  20%  glycerol,  2%  SDS,  50  mM  2- 
mercaptoethanol,  and  0.1%  bromophenol  blue. 
The  protein  samples  were  electrophoresed  on 
4-15%  SDS-PAGE  gradient  gels  (BioRad)  and 
transferred  to  nitrocellulose  membrane.  The 
immunoblots  were  detected  with  anti-schwan- 
nomin  antibody  2781  described  previously 
[Huynh  and  Pulst,  1996].  The  blots  were 
reprobed  with  anti-actin  monoclonal  antibody 
AC40  (Sigma)  for  normalization  of  protein 
loading.  In  RT4  cells,  schwannomin  or  mutant 
schwannomin  was  detected  in  a  similar  manner 
except  that  the  anti-schwannomin  antibody 
WA30  was  used  [Gutmann  et  al.,  1997]. 

RESULTS 

Generation  of  Inducible 
Schwannomin-Expressing  MEF  Cell  Lines 

Inducible  cell  lines  expressing  schwannomin 
protein  isoforms  1  and  2  were  generated  in 
mouse  embryonic  fibroblast  cells  (parental 
MEF/3T3  Tet-Off  lines  obtained  from  Clontech) 
using  retroviral  infection.  The  four  selected  cell 
lines  (4-6,  15-7  for  isoform  1  and  2-19,  23-6  for 
isoform  2)  and  the  parental  cell  line  were 
maintained  in  100-mm  plates  with  2  pg/ml 
doxycycline  in  the  medium.  When  cells  reached 
50%  confluency,  cells  were  washed  thoroughly 
in  37  C  dox-free  media  to  remove  any  trace  of 
doxycycline  and  incubated  again  in  the  media 
with  or  without  doxycycline  for  24  h  to  repress  or 
induce  the  schwannomin  protein  expression. 
The  experiment  was  replicated  for  parental,  4-6, 
15-7,  and  2-19  cell  lines  and  repeated  four  times 
for  line  23-6  (Table  I).  The  repeated  experiments 
were  denoted  with  small  characters,  a,  b,  c,  and 
d  after  cell  line  names.  Based  on  immunoblot¬ 
ting  analysis  (Fig.  1),  several  samples,  such  as 
15-7b,  4-6a,  4-6b,  23-6a,  and  23-6c,  showed  low 
basal  levels  of  schwannomin  expression  and 
strong  induction  of  schwannomin  upon  doxycy¬ 
cline  withdrawal. 

The  purified  RNAs  were  labeled  and  hybri¬ 
dized  for  expression  profiling  using  Affymetrix 
MG_U74Av2  chips.  After  hybridization  and 
image  analysis,  the  image  files  (CEL  files)  were 
imported  into  the  DCHIP  program  [Li  and 
Wong,  2001].  While  calculating  expression 
levels  using  the  PM/MM  difference  model,  the 
program  excluded  the  probe  sets  that  did 
not  match  the  overall  pattern  of  intensities 
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TABLE  I.  Microarray  Experiments  Using  29  Mouse  Embryonic  Fibroblast 
Tissue  Samples  From  Five  Cell  Line  Clones 


Clonal  cell  line 

Array 

name 

P  call 

(%) 

Array  outlier 

(%) 

Single  outlier 

(%) 

Flag 

Parental  cell 

P 

51.5 

0.288 

0.263 

PaD 

50.6 

0.416 

0.260 

PaN 

50.3 

0.625 

0.271 

PhD 

52.1 

0.296 

0.231 

PbN 

48.7 

3.820 

0.903 

* 

MEF  NF2  isoform  1  clone  4-6 

4-6 

52.3 

0.480 

0.266 

4-6  aD 

49.4 

0.400 

0.284 

4-6  aN 

51.5 

0.288 

0.200 

4-6bD 

46.0 

3.692 

1.025 

* 

4-6bN 

53.0 

0.312 

0.234 

MEF  NF2  isoform  1  clone  15-7 

15-7 

49.2 

0.681 

0.417 

15-7aD 

50.2 

0.464 

0.336 

15-7aN 

49.1 

0.496 

0.332 

15-7bD 

51.8 

0.673 

0.244 

15-7bN 

52.2 

0.288 

0.212 

MEF  NF2  isoform  2  clone  2-19 

2-19 

51.7 

0.488 

0.197 

2-19aD 

50.5 

0.248 

0.200 

2-19aN 

51.0 

0.280 

0.294 

2-19bD 

50.2 

0.553 

0.240 

2-19bN 

53.4 

0.320 

0.226 

MEF  NF2  isoform  2  clone  23-6 

23-6 

51.2 

0.208 

0.173 

23-6aD 

52.0 

0.488 

0.339 

23-6aN 

47.7 

1.241 

0.542 

23-6bD 

49.1 

0.416 

0.307 

23-6bN 

46.3 

3.211 

0.952 

• 

23-6cD 

52.8 

0.537 

0.388 

23-6cN 

52.0 

0.569 

0.380 

23-6dD 

52.5 

0.641 

0.334 

23-6dN 

52.1 

0.288 

0.266 

The  cell  lines  used  were  parental  (P),  4-6,  15-7,  2-19,  and  23-6.  Array  namesending  with  N  are  the  samples 
from  24  h  after  the  removal  of  doxycycline,  and  ones  with  D  are  the  control  sample  with  doxycycline  in  the 
medium.  The  small  characters  (a-d)  before  N  or  D  of  the  array  name  represent  the  repeats  of  experiments. 
Samples  without  N  or  D  were  taken  at  0  h  before  removal  of  doxycycline  at  50%  confluence.  P  call  (%) 
represents  percentage  of  the  probes  called  ‘present’  in  the  array  by  DCFIIP  program.  Array  outliers 
represent  the  probe  sets  that  do  not  follow  the  overall  pattern  of  the  intensities  calculated  by  model-based 
expression  tags  in  DCIIIP.  Single  outliers  represent  single  probes  that  do  not  follow  intensity  pattern  of  the 
respective  probe  in  other  sets.  If  the  percentages  of  array  and  single  outliers  are  high,  the  array  is  flagged. 


observed  with  other  chips  or  the  probes  did  not 
match  the  intensity  pattern  of  the  respective 
probe  on  other  chips.  Those  probe  sets  were 
designated  “array  outliers”  or  “single  outliers,” 
respectively  (Table  I).  In  general,  the  outlier 
percentage  is  inversely  correlated  with  the 
image  quality  [Schadt  et  al.,  2000].  Among  the 
29  microarray  experiments,  three  chips,  P2N,  4- 
62D,  and  23-62N,  showed  notably  higher  per¬ 
centages  of  probe/single  outliers  than  others. 
Therefore,  these  three  chips  were  eliminated 
from  subsequent  analyses. 

Effect  of  Schwannomin  Induction 
on  Expression  Profiling 

Schwannomin  expression  was  monitored 
with  immunoblotting  experiments  (Fig.  1). 
Two  samples  each  for  isoform  1  (4-6a,  15-7b) 
and  isoform  2  (23-6a,  and  23-6c)  showed  strong 
induction  and  repression  of  schwannomin 
expression  responsive  to  the  level  of  doxycy¬ 
cline.  Two  pairs  of  expression  profiles,  4-6aD/4- 


6aN  and  15-7bD/15-7bN,  were  analyzed  using 
Student's  £-test  in  the  DCHIP  program  to 
identify  the  genes  regulated  by  schwannomin 
isoform  1  expression.  For  pair-wise  comparison, 
we  first  identified  the  genes  called  'present'  in 
at  least  one  chip  by  the  DCHIP  program.  In 
all  pair-wise  comparison,  approximately  7,000 
genes  were  used  for  further  analysis.  At  signi¬ 
ficance  levels  of  P  <  0.05  and  <  0.01,  one  would 
expect  to  identify  350  and  70  false  positives, 
respectively,  from  a  set  of  7,000  probes.  We 
found  only  one  gene  to  be  significantly  regulated 
in  the  15-7bD/15-7bN  comparisons  and  none  in 
the  4-6aD/4-6aN  comparisons  using  a  signifi¬ 
cance  level  P  <  0.01.  With  a  relaxed  stringency 
(P  <  0.05),  48  and  53  genes  were  filtered  out  as 
changed  in  the  4-6aD/4-6aN  and  the  15-7bD/15- 
7bN  comparisons,  respectively.  In  summary, 
the  number  of  genes  significantly  changed  was 
much  smaller  than  the  expected  number  of  false 
positives,  suggesting  that  expression  levels  of 
very  few  or  no  genes  were  significantly  changed. 
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Fig.  1.  Western  blot  detection  of  schwannomin  expression 
levels  in  MEF/3T3  cell  lines.  Cells  were  incubated  in  the  medium 
with  or  without  doxycycline  for  24  h  followed  by  washing  with 
dox-free  medium.  After  electrophoresis  of  total  cell  lysate, 
schwannomin  was  detected  by  using  antibody  2781  [Huynh 
and  Pu  1st,  1 996 J.  To  control  for  loading,  blots  were  reprobed  with 


an  anti-actin  monoclonal  antibody  AC40  (Sigma).  Detail 
experimental  procedure  is  described  in  the  text.  Parental  cell 
line  samples  show  endogenous  expression  level  of  schwanno¬ 
min.  Two  cell  line  samples,  15-7a  and  23-6d,  show  poor 
suppression  under  the  existence  of  doxycycline  and  marked  with 
asterisk. 


To  detect  genes  differentially  expressed  in 
clonal  cell  lines,  we  also  examined  the  number 
of  regulated  genes  between  the  two  control 
(uninduced)  samples  (4-6aD/15-7bD)  and  be¬ 
tween  the  two  induced  samples  (4-6aN/15-7bN). 
Surprisingly,  these  analyses  found  201  and  159 
differentially  regulated  genes  using  /-test  with 
P<0.05,  respectively.  Among  them,  62  genes 
were  commonly  up-  or  downregulated  in  both 
comparisons.  This  suggests  that  many  genes 
were  differentially  expressed  between  4-6a  and 
15-7b  cell  lines  even  before  the  schwannomin 
isoform  1  was  induced.  The  differential  expres¬ 
sion  profiles  between  different  clonal  cell  lines 
will  be  discussed  below. 

Expression  of  schwannomin  isoform  2  did  not 
result  in  significant  gene  expression  changes. 
According  to  Student’s  /-test  in  the  DCHIP 
program,  none  of  the  probes  was  significantly 
regulated  with  a  value  of  P<0.01  in  compar¬ 
isons  between  the  23-6aN/23-6aD  and  between 
the  23-6cN/23-6cD  chip  sets.  With  a  P-value  < 
0.05,  83  and  18  genes,  respectively,  were  dif¬ 
ferentially  changed,  numbers  smaller  than 
the  expected  false  positives.  Among  them,  five 
genes  were  regulated  in  common,  but  their 
expression  level  changes  were  not  significant 


(<  1.5-fold).  A  similar  number  of  differentially 
regulated  genes  were  found  when  we  compared 
two  replicate  samples  from  the  same  cell  line 
before  and  after  induction.  This  suggests  that 
neither  schwannomin  isoform  caused  signifi¬ 
cant  expression  changes  in  a  large  number  of 
genes  when  overexpressed  in  mouse  embryonic 
fibroblasts. 

Experimental  Design  and  Expression 
Changes  in  RT4  Cells 

The  effect  of  regulated  schwannomin  expres¬ 
sion  was  independently  examined  in  rat  RT4 
schwannoma  cells.  For  these  experiments,  two 
cell  lines,  in  addition  to  the  parental  Tet-On  RT4 
line,  were  used;  one  line  expressed  wild  type 
schwannomin,  the  second  one  schwannomin 
with  the  disease-associated  L64P  amino  acid 
substitution  [Gutmann  et  al.,  2001].  Previous 
studies  have  demonstrated  that  schwannomin 
containing  the  L64P  mutation  is  functionally 
inactive  [Gutmann  et  al.,  2001].  Upon  addition 
of  doxycycline  to  the  media,  schwannomin  was 
consistently  induced  (Fig.  2).  Total  RNAs  were 
purified  from  the  cells  at  0,  6,  and  24  h  after 
induction;  the  experiment  was  replicated  and 
the  samples  designated  by  the  small  letters  a 
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Fig.  2.  Western  blot  analysis  of  protein  expression  in  RT4  cell  line  clones.  Parental  cell  line  showed 
endogenous  schwannomin  expression,  while  cell  lines  transfected  with  schwannomin  and  mutant 
schwannomin  show  significantly  induced  protein  expression  after  adding  doxycycline.  The  proteins  were 
probed  with  WA30  antibody  iGutmann  et  al.,  1997]. 


or  b.  The  expression  levels  of  the  18  samples 
were  monitored  using  the  Affymetrix  RG  U34A 
chips  (Table  II).  In  this  experiment,  MOb  show¬ 
ed  a  high  array/single  outlier  percentage  and 
was  excluded  from  subsequent  analyses. 

The  effect  of  schwannomin  or  mutant  sch¬ 
wannomin  induction  on  gene  expression  pro¬ 
files  was  analyzed  in  the  same  way  as  the  MEF / 
3T3  cell  lines.  The  expression  levels  of  0-h 
samples  from  parental,  schwannomin,  and 
mutant  schwannomin  producing  cell  lines  were 
compared  with  6-  or  24-h  samples  of  the  same 
cell  lines  using  the  /-test.  About  4,500  probe  sets 


were  called  ‘present’  in  more  than  one  chip  in 
each  pair-wise  comparison.  By  the  /-test  with 
P<  0.05,  275  false  positives  would  be  expected. 
The  differentially  expressed  genes  found  by 
the  /-test  between  0-h  samples  and  6-  or  24-h 
samples  did  not  exceed  the  expected  number  of 
false  positives  in  any  comparison.  For  example, 
the  comparison  between  0-  and  24-h  samples 
of  schwannomin  inducible  cell  line  (S0a/S24a 
in  Table  II)  found  34  differentially  expressed 
genes.  This  number  is  even  smaller  than  67 
and  37  genes  found  by  the  replicate  sample 
comparisons  (between  SOa/SOb  and  between 


TABLE  II.  Microarray  Experiments  Using  18  RT4  Cell  Line  Samples 


Clonal  cell  line 

Array 

name 

P  call 

(%) 

Array  outlier 

(%) 

Single  outlier 

(%) 

Flag 

Parental  RT4  cell 

P0a 

46.9 

0.511 

0.100 

POb 

47.3 

0.557 

0.079 

P6a 

49.9 

0.636 

0.073 

P6b 

49.7 

0.443 

0.064 

P24a 

48.3 

0.273 

0.045 

P24b 

49.7 

0.284 

0.071 

Schwannomin  inducible 

SOa 

46.1 

0.443 

0.084 

RT4  cell 

SOb 

46.1 

0.466 

0.104 

S6a 

47.8 

0.784 

0.143 

S6b 

47.7 

0.386 

0.043 

S24a 

48.4 

0.955 

0.130 

S24b 

48.3 

0.341 

0.071 

Mutant  schwannomin 

MOa 

47.7 

0.227 

0.054 

inducible  RT4  cell 

MOb 

44.3 

5.455 

0.870 

* 

M6a 

47.6 

0.841 

0.136 

M6b 

47.6 

1.978 

0.319 

M24a 

47.0 

0.330 

0.073 

M24b 

46.7 

0.580 

0.126 

Array  names  starting  with  P,  S,  and  M  represent  parental,  schwannomin-,  and  mutant  schwannomin- 
ex pressing  cell  lines,  respectively.  The  numbers  following  the  letter  represent  the  sampling  time  after  the 
induction  of  protein  expression,  and  a  and  b  the  replication  of  the  experiment.  The  remainder  of  the  terms 
are  the  same  as  in  Table  I. 
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S24a/S24b),  which  indicates  that  schwannomin 
expression  had  little  effect  on  gene  expression 
profiles  in  RT4  cells. 

Hierarchical  Clustering  of  the  Data 

Instead  of  comparing  differentially  expres¬ 
sed  genes  by  /-test,  the  overall  pattern  of  the 
expression  profiles  was  investigated  with  clus¬ 
ter  analysis  in  MEF  cells.  Among  12,300  probe 
sets  on  the  Affymetrix  MG_U74Av2  chip,  7,185 
probe  sets  were  chosen  after  eliminating  the 
probe  sets  called  ‘absent’  in  more  than  80%  chips 
by  DCHIP  among  the  26  chip  data.  The  expres¬ 
sion  levels  of  those  probe  sets  were  analyzed  by 
an  unsupervised  cluster  analysis.  The  expres¬ 
sion  profiles  formed  five  major  clusters.  Each 
cluster  contained  the  samples  originating  from 
the  same  clonal  cell  lines  (Fig.  3a).  This  in¬ 
dicated  that  samples  expressing  schwannomin 
did  not  cluster  together,  but  instead  that  clonal 
selection  during  the  generation  of  stably  expres¬ 
sing  cell  lines  had  resulted  in  distinct  expres¬ 
sion  profiles. 

The  cluster  analysis  was  re-examined,  this 
time  focusing  on  probes  with  high  expression 
levels  that  are  known  to  be  more  reliable 
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[Tu  et  al.,  2002].  We  chose  958  probe  sets 
with  intensities  over  500  (maximum  inten¬ 
sity  ~  6400)  and  presence  in  more  than  80%  of 
chips.  The  hierarchical  clustering  revealed  a 
similar  cluster  tree  pattern  to  the  one  with  7, 185 
probe  sets,  each  demonstrating  five  clusters 
with  samples  from  different  cell  lines.  Minor 
exceptions  were  identified:  the  sample  4-6aD 
clustered  with  the  15-7  cell  line  samples  and 
the  sample  P  clustered  with  the  sample  2-19 
(Fig.  3b). 

Cluster  analysis  with  the  expression  profiles 
of  RT4  cell  lines  showed  similar  results.  Among 
the  8,700  probes  on  the  RG  34A  chip,  4,689 
probe  sets  were  called  ‘present’  in  more  than 
20%  of  the  chips.  The  cluster  analysis  with  the 
expression  profiles  of  the  probe  sets  showed  that 
17  samples  formed  three  major  clusters  with 
each  cluster  containing  the  samples  from  the 
same  clonal  cell  line  (Fig.  4).  The  time  points 
after  the  protein  induction  had  few  effects  on 
global  gene  expression  levels,  although  the  6-  or 
24-h  samples  expressed  significant  levels  of 
wild  type  or  mutant  schwannomin  expression 
(Fig.  2).  This  ascertained  that  the  variations 
incorporated  during  selection  process  were  more 
significant  than  the  ones  induced  by  schwanno¬ 
min  protein  overexpression. 

Differentially  Expressed  Genes 
in  Clonal  Cell  Lines 

To  identify  an  even  smaller  subset  of  differ¬ 
entially  expressed  genes  in  different  clonal  cell 
lines,  Kruskal-Wallis  rank  test  was  performed 
with  log- transformed  gene  expression  ratios. 


Fig.  3.  Unsupervised  hierarchical  clustering  of  expression 
profiles  of  26  mouse  embryo  fibroblast  cell  line  samples  listed 
in  Table  I.  The  cluster  analysis  was  performed  with  (a)  7,185 
probe  sets  and  (b)  958  probe  sets  selected  based  on  the 
percentages  of  present  calls  and/or  expression  levels  as  described 
in  the  text.  The  length  of  vertical  bars  indicates  the  degree  of 
difference  in  gene  expression  levels  between  samples.  The 
samples  from  the  same  cell  line  clone  generally  cluster  together 
irrespective  of  the  induction  of  schwannomin  protein  expression. 


Fig.  4.  Unsupervised  hierarchical  clustering  of  gene  expression 
data  of  1  7  RT4  cell  lines  in  T able  II.  The  clustered  tree  shows  that 
different  cell  line  clone  is  the  major  variable  that  makes 
differences  in  the  expression  levels.  Among  the  cell  line  clones, 
parental  and  schwannomin  transfected  cell  lines  show  similar 
expression  levels  and  missense  mutant  schwannomin  transfected 
cell  line  shows  the  most  distinct  expression  profile. 
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The  expression  levels  of  the  samples  from  the 
same  clonal  cell  line  were  grouped  together, 
generating  five  and  three  groups  for  experi¬ 
ments  with  MEF/3T3  and  RT4  cell  lines, 
respectively.  The  genes  differentially  expressed 
in  different  groups  were  searched  for  with 
Kruskal-Wallis  rank  test.  The  analyses  with 
MEF/3T3  and  RT4  cell  lines  selected  320  and 
292  genes,  respectively,  as  “highly  expressed” 
(intensities  over  500  and  ‘present’  in  more  than 
80%  of  the  chips)  and  as  “significantly  varied” 
(P<0.01  in  Kruskal-Wallis  rank  test)  in  dif¬ 
ferent  clonal  cell  lines.  A  total  of  28  genes 
variable  in  both  mouse  and  rat  cell  lines  were 
identified  in  this  fashion.  Of  these  transcripts, 
13  genes  were  identical  for  both  cell  types  and 
15  were  highly  related.  Of  the  28  genes,  a 
significant  number  represented  cytoskeletal 
related  genes  (8  out  of  28)  including  actin, 
annexin  A1  (lipocortin-1),  fibronectin,  integrin, 
osteonectin,  and  tubulin.  The  remaining  genes 
included  cytochrome-c  oxidase  subunits,  pro- 


teasomal  proteins, and  ribosomal  proteins.  Most 
of  the  28  genes  varied  less  than  2-fold  among  the 
cell  lines  with  the  highest  and  lowest  measure¬ 
ments.  However,  a  few  genes,  such  as  proteaso- 
mal  subunit  a  and  protein  phosphatase  4,  were 
differentially  expressed  up  to  threefold  among 
the  RT4  cell  line  clones.  The  list  of  28  genes  and 
their  expression  levels  are  shown  in  Figure  5. 

DISCUSSION 

Subtle  Expression  Level  Changes 
by  Schwannomin  Inductions 

We  performed  microarray  experiments  with 
two  stably  transfected  cell  lines,  MEF/3T3  and 
RT4,  expressing  schwannomin  under  the  con¬ 
trol  of  an  inducible  promoter.  Our  experiments 
could  not  identify  significantly  regulated  genes 
induced  by  schwannomin  overexpression.  This 
was  unexpected,  considering  that  schwanno¬ 
min  is  a  tumor  suppressor  protein  involved  in 
signal  transduction  pathways  mediated  by  Rho/ 
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Fig.  5.  The  list  of  the  genes  varied  significantly  in  different  MEF/ 
3T3  and  RT4  cell  line  clones.  The  genes  were  chosen  using 
Kruskal-Wallis  rank  test  as  described  in  the  text.  Bright  red  and 
green  colors  represent  the  expression  level  more  than  50%  above 
and  below  the  average  level  in  all  samples,  respectively.  The 
darker  the  color,  the  lesser  the  change  in  expression  level.  The  left 
clustering  tree  is  the  expression  profile  of  MEF/3T3  cell  line 
samples,  while  the  right  one  is  of  RT4  cells.  The  samples  from  the 


same  clonal  cell  lines  were  clustered  together  except  the  sample 
2-19.  Among  the  genes  presented  in  this  figure,  13  are  the 
functionally  identical  genes  in  both  cell  lines  and  the  other  15 
with  asterisk  at  the  end  of  gene  name  are  the  genes  with  similar 
function  according  to  the  mouse  and  rat  gene  annotations.  The 
first  term  in  the  parenthesis  is  the  gene  name  corresponding  to 
MEF/3T3  cell  and  the  second  to  RT4  cell. 
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Rac  and  STAT3/5  IPelton  et  al.,  1998;  Scoles 
et  al.,  2002].  It  is  possible  that  changes  in 
expression  levels  were  so  subtle  that  they 
were  not  detected  in  chip  hybridizations  or  that 
the  relevant  genes  were  not  represented  on  the 
chips  used  in  our  experiments.  In  addition, 
gene  expression  changes  may  have  been  mis¬ 
sed,  because  they  were  asynchronous  or  occur¬ 
red  outside  the  time  frames  that  we  examined. 
It  is  also  possible  that  the  major  effects  of 
schwannomin  overexpression  result  in  direct 
protein  modification  without  significantly  alter¬ 
ing  steady-state  gene  expression  profiles. 

Another  possibility  is  that  transcriptional 
changes  induced  by  schwannomin  are  only 
noticeable  under  certain  experimental  condi¬ 
tions  in  vitro.  For  example,  we  previously 
showed  that  schwannomin  regulated  IGF-I 
induced  STATS  phosphorylation  after  serum- 
starvation  [Scoles  et  al.,  2002].  Phosphorylated- 
STAT3  binds  to  DNA  and  alters  expression 
levels  of  many  genes.  The  experimental  condi¬ 
tion  used  in  this  study  may  not  have  been 
appropriate  to  detect  those  physiological 
changes.  Finally,  whereas  the  loss  of  schwan¬ 
nomin  protein  causes  significant  physiological 
effects  such  as  the  formation  of  tumors  in  the 
peripheral  nervous  system,  the  induction  of 
schwannomin  protein  may  not  cause  any 
severe  physiological  effects  in  cells  maintained 
in  culture.  Indeed,  cDNA  microarray  data  with 
human  schwannomas  showed  dramatic  tran¬ 
scriptional  changes,  such  as  upregulation  of 
osteonectin  and  RhoB  GTPase  as  well  as  down- 
regulation  of  LUCA-5  and  CDK2  [Lasak  et  al., 
2002;  Welling  et  al.,  2002].  Expression  level 
changes  in  human  schwannomas  in  those 
studies  require  further  verification.  Small  sam¬ 
ple  size,  tissue  heterogeneity,  and  inter-indivi¬ 
dual  variations  among  human  patients  may 
result  in  a  substantial  rate  of  false  positives 
[Bakay  et  al.,  2002]. 

Microarray  experiments  have  not  always 
been  able  to  identify  gene  expression  changes 
despite  the  presence  of  significant  phenotypes. 
Detailed  microarray  experiments  with  mouse 
brains  lacking  the  MECP2  gene  found  no  signi¬ 
ficant  transcriptional  changes  [Tudor  et  al., 
2002].  The  result  was  even  more  remarkable,  as 
MECP2  is  known  as  a  transcriptional  repres¬ 
sor  and  its  mutation  causes  a  severe  cognitive 
phenotype  in  humans,  Rett  syndrome.  Our 
microarray  experiments  strongly  suggest  that 
schwannomin  overexpression  in  two  different 


cell  types  caused  few  transcriptional  changes 
as  well. 

Heterogeneity  Introduced  by 
Clonal  Selection 

Transfection  of  a  gene  of  interest  into  cultur¬ 
ed  cells  is  one  of  the  most  commonly  used 
methods  to  study  the  cellular  function  of  the 
respective  protein.  The  high-throughput  scre¬ 
ening  ability  of  microarray  experiments  has 
provided  a  powerful  tool  to  characterize  protein 
function  in  transfected  cell  lines.  Various  ex¬ 
perimental  designs  have  been  used,  such  as 
transient  transfection  or  stable  transfection 
with  inducible  or  constitutive  promoters.  How¬ 
ever,  the  background  variability  introduced  by 
transfection  or  the  selection  process  has  not 
been  systematically  studied. 

With  microarray  experiments  of  stably  trans¬ 
fected  cell  lines,  we  observed  significant  altera¬ 
tions  of  gene  expression  between  different 
clonal  cell  lines.  The  unsupervised  cluster 
analyses  with  26  MEF/3T3  and  1 7  RT4  cell  line 
samples  indicated  that  clonal  selection  was  the 
major  variable  that  changed  gene  expression 
levels  independent  of  expression  of  exogenous 
schwannomin.  Although  technical  variations 
could  confound  the  results,  they  would  not 
contribute  much  to  the  observed  variations, 
because  in  each  experiment,  RNA  purifications 
were  carefully  performed  under  the  same  con¬ 
ditions  and  all  chip  experiments  were  per¬ 
formed  on  the  same  day  in  random  orders.  In 
addition,  replicate  samples  from  the  same  cell 
line  and  treatment  condition  (indicated  by  small 
letters  in  Figs.  3  and  4)  clustered  together 
indicating  a  high  degree  of  reproducibility  of 
the  data. 

We  used  Kruskal -Wallis  rank  test  to  identify 
a  small  subset  of  genes  that  were  highly 
variably  expressed  across  all  different  clones 
including  mouse  and  rat  cell  lines.  This  small 
set  of  28  genes  resulted  in  a  tree  that  was  highly 
similar  to  the  ones  generated  with  a  larger 
number  of  genes.  Inspection  of  the  expression 
patterns  of  single  genes  indicated  that  expres¬ 
sion  was  highly  different  between  different  cell 
line  clones.  For  example,  annexin  Al  expression 
was  high  in  23-6  cells,  but  low  in  mouse  parental 
cells  and  in  15-7  cells.  The  inverse  was  observed 
for  peroxiredoxin  that  showed  higher  expres¬ 
sion  in  all  15-7  clones  compared  with  23-6  clones 
independent  of  schwannomin  expression. 
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The  expression  of  many  cytoskeletal  proteins 
varied  in  expression  between  the  different 
clones.  It  should  be  noted  that  some  of  these 
proteins,  such  as  actin  and  tubulin,  are  being 
used  as  reference  proteins  to  normalize  the 
gene  or  protein  expression  levels  in  Northern  or 
Western  blots.  Our  results  strongly  suggest  that 
the  expression  levels  of  these  genes  may  vary  in 
different  clonally  derived  cell  lines,  although 
the  level  changes  were  less  than  40%  between 
the  clones. 

The  changes  in  background  expression  may 
not  be  relevant  when  examining  the  effects 
of  exogenously  expressed  genes  that  have  a 
strong  effect  on  gene  expression.  For  example,  a 
number  of  regulated  genes  as  a  result  of  over¬ 
expression  of  BRCA1  [Welcsh  et  al.,  2002)  or 
mutant  huntingtin  [Sipione  et  al.,  2002]  were 
identified  despite  background  variability  of 
stably  transfected  cells.  On  the  other  hand, 
clonal  background  changes  may  greatly  influ¬ 
ence  results,  when  exogenously  expressed 
transgenes  have  weak  effects,  because  expres¬ 
sion  levels  for  some  genes  may  greatly  differ 
between  cell  clones  that  appear  phenotypically 
highly  similar. 

The  reason  why  substantial  expression  level 
changes  are  introduced  during  transfection  or 
selection  process  requires  further  study,  but 
may  relate  to  the  chromosomal  location  of 
plasmid  integration  or  the  selection  of  particu¬ 
lar  cells  in  the  transfected  pool  with  pre-existing 
subtle  differences  in  background  gene  expres¬ 
sion.  This  background  variability  should  be  con¬ 
sidered  during  experimental  design  and  data 
interpretation.  In  the  mean  time,  our  results 
suggest  that  microarray  experiments  with 
stably  transfected  genes  should  be  replicated 
with  multiple  cell  line  clones  and  that  the  direct 
comparison  of  expression  data  from  different 
clones  should  be  interpreted  with  caution. 
Alternatively,  transient  transfection  of  the  gene 
can  reduce  the  background  variability  if  the 
transfection  method  is  highly  efficient  or  if  it  is 
combined  with  a  selection  method. 
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